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ABSTRACT 

We investigate the effect of inverse-Compton scattering of ffares of soft 
radiation in different geometries of a liot, Comptonizing region and a colder 
accretion disk around a solar-mass black hole. The photon-energy dependent 
light curves, their Fourier transforms, power spectra and Fourier-period 
dependent time lags of hard photons with respect to softer photons are 
discussed. On the basis of a comparison with existing data we find arguments 
against Comptonization of external soft radiation as well as Comptonization 
in a homogeneous medium as dominant mechanisms for the rapid aperiodic 
variability in Galactic black-hole candidates. Possible further observational 
tests for the influence of Comptonization on the rapid aperiodic variability of 
Galactic black-hole candidates are suggested. 

Subject headings: X-rays: stars — accretion, accretion disks — black hole 
physics — radiative transfer — radiation mechanisms: thermal 
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1. Introduction 

The rapid aperiodic and quasiperiodic variability of the X-ray signals from Galactic 



black-hole candidates (GBHCs) and low-mass X-ray binaries (for a review see 



van aer 



[Klis 19951) contains valuable information about the source of high-energy emission in these 



objects. Its typical timescales, typical repetition frequency, power spectrum density (PSD), 
autocorrelation function, time lags between different energy bands etc. give hints toward 
important parameters such as the extent of the emitting region, the dominant microscopic 
timescales, the dominant emission mechanism and the geometry of the source. 

With the new generation of X-ray telescopes such as the PGA and the HEXTE on 
board RXTE it is now possible to measure the above properties of the rapid variability of 
X-ray binaries in great detail with high timing resolution and good spectral resolution. As 
we will see in this paper, the photon-energy dependence of the rapid aperiodic variability 
can provide important diagnostics of the nature of the X-ray emitting regions. 

Early measurements of the Fourier-frequency dependence of time lags between the 
signals in different X-ray energy channels from GBHGs ( [Miyamoto et al. 1988| , [Miyamoto 



et al. 19"93D have been interpreted as evidence that Gomptonization of soft photons in a 



hot, uniform plasma could not be the dominating mechanism for the production of hard 
X-rays. This conslusion, however, is strongly geometry-dependent and does not hold for a 
very extended, inhomogeneous Gomptonizing region, as was recently found by Kazanas et 
al. (1997). In that and in two subsequent papers ( |Hua et al. 19974 [Hua et al. 1997b|) the 



effect of radial density gradients in the hot Gomptonizing regions in GBHGs on the power 
density spectra and the phase and time lags between different energy bands were discussed. 
From the comparison of the measured Fourier-frequency dependent time lags between two 
energy channels of the signal from Gyg X-1, they deduced that a radial density dependence 
n(r) oc in the Gomptonizing hot plasma is appropriate to account for the observed hard 
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time lags and the hard X-ray spectrum at the same time. 

However, their work was restricted to a geometry with a central soft photon source 
surrounded by a spherical corona of hot Comptonizing plasma. Furthermore, the density 
was the only parameter allowed to vary radially. 

In this paper, we extend the investigation of the effects of Comptonization the rapid 
variability in X-ray binaries, discussing two fundamentally different source geometries and 
radial gradients of other physical parameters, such as the electron temperature. This 
is done primarily with Monte-Carlo simulations of instantaneous flares of soft radiation 
(described as a 5-function in time), being Comptonized by the X-ray emitting region. 
The generalization to the case of multiple, randomly distributed flares (shot-noise) is 
straightforward. The light curves resulting from the Monte-Carlo simulations are then 
Fourier transformed, and the PSD as well as the Fourier-frequency dependent phase and 
time lags between different energy bands are calculated. 

In Section 2, we describe our Monte-Carlo simulations leading to the energy-dependent 
light curves for the different geometrical situations, together with analytical approximations 
for these light curves. Their Fourier transforms and the corresponding time lags between 
different energy bands will be discussed in Section 3. In Section 4, we investigate specific 
differences between the PSDs and hard time lags resulting from different geometries and 
compare the predictions of both scenarios to the power spectra and hard time lags observed 
for some GBHCs. We will suggest further observational tests which could either confirm a 
preferred geometry, or rule out to the Comptonization scenario in general. 
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2. Simulated light curves 

For our time-dependent Monte-Carlo simulations on Compton scattering in a hot 
plasma we use an improved version of the Comptonization Code developed by Canfield et al. 
(1987), see also Liang (1993), Wandel & Liang (1991). It can handle spherically symmetric 
or slab geometry of the Comptonizing region with thermal soft photon sources located in the 
center of the Comptonizing cloud, distributed throughout the cloud, and/or located outside 
the Comptonizing region. All parameters (such as density, electron temperature, radiation 
temperature of an internal soft photon field) may have arbitrary radial dependencies within 
the cloud. 

Practically, we split the Comptonizing region into 20 radial zones, each zone having 
the same (radial) Thomson depth, within which the parameters are set to a constant value 
corresponding to an average of the respective radial dependencies over the zone. For several 
of our simulations, we run test simulations with a finer radial grid (50 - 100 zones) in order 
to check that the step function dependency of the parameters in the Comptonizing region 
does not introduce artificial structures in the light curves. No significant difference to our 
simpler (and faster) 20 zone simulations has been found in any of these tests. 

In this paper, we focus on two basic situations with respect to the location of the flaring 
soft photon source: a source located at the center of the cloud and a source located outside 
the cloud. The former situation is representative of an accretion-disk-corona model ([Liang 



[and Price 1977| , [Bisnovatyi-Kogan and Blinnikov 1977| , [Haardt and Maraschi 1993|) where 



photons resulting from a flare in the cold, optically thick disk are Compton-upscattered 
by the hot corona to produce the hard X-ray spectrum. The latter case represents a 
hot-inner-disk model ( [Shapiro et al. 1976| ) or an advection-dominated accretion flow 
( [Narayan and yi 1994[ , [Chen et al. 1995| ), assuming that a soft photon flare occurs in the 



cold, optically thick outer disk and its radiation impinges on the hot inner disk (torus) 
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producing the hard X-rays via Compton upscattering. 

For all simulations shown in this paper, the soft photon spectrum is assumed to be a 
thermal blackbody of kTr — 0.2 keV. 

For each of the light curves resulting from central soft photon injection, a total of 9 • 10^ 
photons, being injected instantaneously at t — 0, are simulated. The photons leaving the 
system are collected in 5 energy channels and in equal time steps appropriate to reveal the 
detailed structure of the light curves (typically At ~ 10~^ i?/c, where R is the radius of the 
Comptonizing region). The light curves resulting from external soft photon injection turn 
out to be mathematically more complicated and require a higher numerical signal-to-noise 
ratio in order to yield reliable Fourier transforms. Therefore, the simulations for this 
geometry have been carried out with 25 million photons each. 

For both geometries, we investigate the effect of density and temperature gradients 
within the hot Comptonizing region. Fig. 1 shows the energy-dependent light curves for 
the case of a central photon source inside a homogeneous cloud. We did further simulations 
for Comptonizing clouds with density profiles n{r) oc for p = 1 and 3/2 and for clouds 
with temperature profiles of the same form. Generally, the light curves for a central source 
can be parametrized by a Gamma distribution. 



as suggested by Hua et al. (1997b), where 




if X < 

1 if X > 



is the Heaviside function and ^4 is a normalization constant. Inspecting the simulated light 
curves for the different cases, we find that both a and (3 are energy- dependent in both 
the case of a homogeneous cloud as well as a cloud with parameters spatially varying. 
This is in contradiction to the finding of Hua et al. (1997b) who found a qualitative 




(1) 
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difference between the light curves in these two cases. According to our simulations, 
the difference between these two situations is purely quantitative. The choice of a much 
smaller homogeneous Comptonizing region compared to the case of an extended cloud with 
n{r) oc r"^ leads to the fact that the differences in the short-term behavior {t <^ 1/P) for 
different energy channels only affects very high Fourier frequencies for which no phase and 
time lags are measured and which have not been resolved in the simulations of Hua et al. 
(1997b). However, in the situation investigated there, this effect only affects very high 
Fourier frequencies for which no phase and time lags are measured anyway. 

Both the early-time slope a, describing the light curve immediately after the onset of 
the flare (accounting for the hght travel time) and the cut-off time P (which is of the order 
of the photon escape time tgsc ~ ttR/c for tt >,1) tend to increase with photon energy. 
The assumption of a density gradient (conserving the total Thomson depth tt and extent R 
of the region) leads to a considerable reduction of a and a moderate reduction of p. These 
changes are much more pronounced for the high photon energies than for the lower energy 
channels, implying that the energy dependence of a and f3 is reduced due to the spatial 
gradient. The assumption of a temperature gradient has a similar effect on the light curve 
parameters than a density gradient. In both cases, scattering in the outer regions of the 
cloud becomes less efficient, leading to a steeper light curve. 

The values of a and P yielding a good representation of the simulated hght curves 
for central soft photon injection are listed in Table 1. However, we point out that in the 
case of radial parameter gradients the light curves of the higher two energy channels show 
significant deviations from a power-law at the onset of the flare of Comptonized radiation. 
This issue will be discussed further in the next section. 

For the case of radiation impinging from outside onto the Comptonizing region, we 
assume that the source of soft photons is located close to the outer boundary of the hot 



8 



plasma region, implying that the radiation enters the Comptonizing region at angle cosines 
with respect to the radial direction which are randomly distributed between -1 and 0. In 
this case, the light curve can generally be parametrized by the sum of a power-law resulting 
from single scattering over the light crossing time through the source, 2R/c plus a multiple 
scattering component, described by a hght curve of the form of Eq. (1), i. e. 



where A and B are normalization constants. Figs. 2 and 3 illustrate the energy-dependent 
light curves for scattering of an external soft radiation flare in a homogeneous and an 
inhomogeneous cloud, respectively, of the same total radial extent and Thomson depth. The 
resulting Comptonization spectra of these two situations are similar to the ones obtained 
from the situations assumed in the case of central photon injection. 

For a homogeneous cloud, clearly the first term in Eq. (3) is dominant, i. e. the 
parametrization is basically the same as in the case of a central soft photon source, except 
in the two highest photon energy channels, where the light curves at t <^ R/c are given by 
a broken power-law with a flat early-time slope. With a more pronounced radial density 
gradient, the second term in Eq. (3) becomes more important and dominates the shape of 
all light curves until t ^ R/c for density profiles with p >l, where n{r) oc r""^. However, 
although the parametrization of the light curves for a uniform cloud scattering radiation 
impinging from the outside is formally the same as for scattering of radiation coming from a 
central soft photon source, the energy-dependence of the spectral index a of the light curves 
for t <C R/c is significantly weaker in the case of the external source. This becomes obvious 
when comparing Tables 1 and 2, where the parameters characterizing the light curves for 
the two geometrically different situations are listed, and from inspection of Fig. 10, where 
the energy-dependence of the early-time slope of the light curves is plotted. 

It should be noted that the light curves resulting from our numerical simulations differ 




(3) 
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significantly from the analytical results of Lightman and Rybicki (1979b) for the problem 
of inverse-Compton reflection off a hot plasma. This descrepancy is due to the fact that 
Lightman and Rybicki used the photon escape probability corresponding to the case of a 
semi-infinite, plane-parallel medium with infinite Thomson depth. Lightman and Rybicki 
(1997a) have shown that the inverse-Compton spectrum in that case is basically given by a 
power-law of photon index 1, which is much harder than the observed hard X-ray spectrum 
of Galactic black- hole candidates (typically F ~ 1.5 - 2). Therefore we have to consider a 
medium of finite Thomson depth, tt <^ 10. 

3. Fourier transforms, power spectra and time lags 

We performed a numerical Fourier Transform for each energy band in our simulated 
light curves and calculated the resulting power spectrum densities. Figs. 4 and 5 show 
the PSD corresponding to the simulations for external soft photon injection for the cases 
p — and p — 1, respectively, illustrated in Figs. 2 and 3. Knowing the Fourier transform 
Fkijjj) = TkC^'^'' of the light curve in energy channel k at Fourier frequency ui = 2%!/, the 
phase difference A4>ki between two energy channels k and / can be computed. The resulting 
time lag between the signals in different energy channels is then given by Atki — A(f)]^i/uj. 

The time lags as a function of Fourier period P = l/i> resulting from our simulations 
for central soft photon injection into a homogeneous and a p — 1 cloud, respectively, are 
displayed in Figs. 6 and 7. Figs. 8 and 9 show the respective time lag curves for the case 
of external soft photon injection. In order to suppress numerical noise contained in the 
simulated light curves, we use smooth analytical curves following the light curves closely 
for the computation of the hard time lags. In Figs. 4 and 5, the PSDs from our analytical 
representations (thin curves) are shown in comparison to the PSDs of the simulated light 
curves (thick curves). The figures demonstrate that the simulated signal of the highest 
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energy channel (with the worst photon statistics) is strongly dominated by numerical noise 
at high frequencies. We therefore do not consider our results concerning the highest energy 
channel (31 - 70 keV) in the case of external photon injection as fully reliable. 

Generally, for both the central and the external soft photon source, the time lags are 
increasing with Fourier period for frequencies uj >, j3^^ ~ c/ {ttR) in the case of a central 
soft photon source and for frequencies uj '>^ujcr = c/{2R) for an external soft photon source. 
For lower frequencies (longer periods) the respective time lags turn into a flat curve where 

For an external soft photon source, no significant deviation from linearity is found, 
independent of parameter gradients or photon energy, except at the lowest photon energies 
which are close to the energy of the injected soft photons. However, in the case of a central 
source the gradual rise of the hght curves at the onset of the flare, systematically deviating 
from a power law, causes a flattening of the time lag dependence on the Fourier period with 
increasing photon energies of the respective channels. In Fig. 6 (central source, p = 0), the 
time lag between the lowest two energy channels depends on the Fourier period to a power 
greater than unity, which we attribute to the fact that the early-time slopes of the two light 
curves in these energy channels are nearly equal (see below). 

The effect of the flattening of the hard-time-lag curve with increasing photon energy 
becomes more pronounced with a steeper density gradient. For a density n(r) cx r~^/^, 
especially the hard lag between the lower energy channels deviates from a simple power-law; 
for the higher channels the curves show power-law behaviour with a slope smaller than 
unity. In this case, the long-period behavior is no longer flat, but the hard time lag begins 
to decline with Fourier period for frequencies lu <, (3~^. 

The power spectra and hard time lags resulting from our simulations show a 
significantly different energy dependence for the two different geometries. In the case of the 
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central soft photon source and a density or temperature gradient, the PSD is flat for low 
photon energies, while for higher photon energies it becomes gradually steeper to higher 
Fourier frequencies. For the p = 1 case, this can be approximated by a broken power-law 
with break frequency cobr ~ P~^- For the solid sphere, the PSD can generally well be 
described by a broken power-law with flat low-frequency slope and the high-frequency slope 
increasing with photon energy. 

For an external soft photon source the PSD can roughly be represented by a 
broken power-law for all energy channels with break frequency Ubr ~ c/ (2i?) and a fiat 
low-frequency branch. The power-law slope depends only very weakly on photon energy 
and is always close to 2. 

Both the analytical representations (1) and (3) of the light curves in the two geometrical 
cases considered here can be Fourier-transformed analytically. The Fourier transforms are 



p. J,.\ _ A r(Qi)/3 iaarctan(/?a;) 

^ " (1 + [/3^]T/^ 



(4) 



for the central soft photon source and 



Fe^tioo) = ^l''^Ca/2 e^"^^*^"(^") + e-'^l^- M f 1; « + 1; i—\ , (5) 

(1 + ^ \ LOcrJ 

where M{a\h\z) is Rummer's function, for the external source. In the limiting cases 
ijj <S P~'^,u!ar and cu 3> (5~^,ujcr-i respectively, these Fourier transforms are 

FinM ^{ (6) 
[AT{a)uj-"e'"''/'^ for a; > 

and 

f A Via) e^"^'^ + e^'^/'^- for uj ^ ujcr 

FextH^l " (7) 

(at (a) uj-"^ e^"'^/^ + B T{k) uj-'^ e^'^'^/^ ^j. ^ ^ 

For the central soft photon source, this yields a Fourier-frequency independent power 
spectrum in the low-frequency limit, while for high frequencies PSDint oc a;~^". In general. 
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the hard time lags between two energy bands in the high- and low-frequency limits are 

I 



r A{al3) for u < 
^tint = { (8) 
AaP/4 fora;>/?-^ 



as found by Hua et al. (1997b). If the values of a for two energy bands are exactly equal, 
then it is easy to show that the phase difference in the high-frequency limit is Acp^i oc a;~^ 
and Atki^aAppy{A7r^/3k/3i). 

From our simulations, we find that the Fourier transforms of the light curves in the 
case of an external photon source are dominated by the second term in Eqs. 5 and 7, 
respectively, in both limits cu <^ cUcr and cu ^ cUcr, if the scattering medium has a strong 
radial density gradient (p >,0.5). The high-frequency hmit of the Fourier transform (5) 
is completely determined by the shape of the light curves at t <^ R/c and can always be 
obtained by simply Fourier transforming a function for the form f{t) — t^~^ describing 
the light curve for t <^ R/c, without accounting for the cutoff at i fa 2R/c or t — P, 
respectively. Therefore, for u ^ Ucr we have PSD oc cu"^'^ where = a for a weak density 
gradient (p^O.5) and low photon energies and /i = k else. 

However, as mentioned in the previous section, this relation between the early-time 
slope of the light curve and the high-frequency branch of the PSD does not apply to the 
case of the free-fall density gradient p — 3/2 and internal soft photon injection, where the 
light curves at the onset of the flare significantly deviate from a power-law behavior, which 
leads to a deviation of the PSD from a simple power-law at high Fourier frequencies. 

For the hard time lags we find in the case of external soft photon injection 

{— for a; ^ LUcr 

(9) 
TT^ P A/I for CU » CUcr- 

The qualitative similarity of the asymptotic expressions for the Fourier transforms for 
both injection geometries is a consequence of the fact that in both cases the rise of the 
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Compton-scattered flare is described by a power-law, which results in a power-law PSD at 
high frequencies, while for both geometries there is a fast decay after one light-crossing time 
tc, which always results in a constant PSD at frequencies lu ^ t~^. 

Assuming that the real signal we measure from an X-ray binary is composed of multiple 
shots randomly distributed in time and with a certain distribution of shot amplitudes An, 
i. e. 

oo 

f{t)= Yl AnUt-Q (10) 

n=— oo 

where /o is one of the single-shot light curves (1) or (3), yields a Fourier transform 

FH=Fo(a;) £ A^e"^-*" (11) 

n=— oo 

where Fq is the respective single-shot Fourier transform (4) or (5). From Eq. (11) it 
is obvious that the time lags between different energy bands remain unaffected by the 
summation of multiple shots, and in the PSD only near the average shot repetition 
frequency some additional structure is introduced. Therefore the asymptotic expressions 
(6) and (7) yield a very good description of the expected high-frequency branch of the PSD 
from a model of Comptonized random shots of soft photons in the two geometries discussed 
here, and Eqs. (8) and (9) remain valid without modification. 



4. Comparison to observations 

The analysis of the frequency- dependent light curves and their Fourier transforms 
enables us to extract important information from observed or observable properties of the 
rapid aperiodic variability of Galactic black-hole candidates. Observable quantities on 
which we will focus now arc: a) the high-frequency power-law index of the PSD and its 
photon-energy dependence, b) the turnover frequency of the PSD, c) the photon-energy 
dependence of hard time lags with respect to the soft X-rays, d) the turnover frequency in 
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the hard-time-lag versus Fourier-Period relation. 

As we have seen in the previous section, the power-law index of the PSD at high 
Fourier frequencies is dominated by the short-term behaviour of individual flares and is 
described by the parameter a or k of Eqs. (1) or (3), respectively. This dependence can be 
written as PSD oc uj^"^^ where /i is the dominant short-time light curve slope (a or k). In 
Fig. 10, we have plotted this slope as a function of the mean photon energy in the energy 
bins in which the light curves were sampled, for density gradients of p = 0, 1, and 3/2, 
respectively, for both central and external soft photon injection. As we mentioned earlier, 
in the case of central soft photon injection and a density gradient of p 1, a significant 
photon energy dependence results, while for an external soft photon source the curves are 
consistent with a constant slope (recalling that we do not attribute high significance to our 
results concerning the highest energy channel for external soft photon injection). 

For the external soft photon source, we find on the basis of our simulations ^ext ~ 1, 
irrespective of radial parameter gradients. This yields a PSD declining as PSD oc cu"^, 
which is much steeper than the observed power spectra which are generally fiat below a 
frequency of z/ ~ 0.1 Hz and turn into a power-law PSD oc v^^ . The slope F has been 
measured for Cyg X-1 to be F ?5i 1 ([Brinkman et al. 1974| , [Nolan et al. 19811) , while for 



GX 339-4 Motch et al. (1983) found F = 1.6 ± 0.2, and Maejimo et al. (1984) evaluated 
F = 1.7. The universality of the external-source PSD slope a 2 therefore are marginally 
consistent with the value found for GX 339-4, but seems to argue strongly against this 
geometry in the case of Cyg X-1. 

For the case of the central source and a density gradient p = 1, our values of a, 
which are in agreement with the ones found by Hua et al. (1997 b), are consistent with 
a PSD oc cu"^ for the energy channel 2.2 — 4.6 keV. For lower energies, our simulations 
predict a slightly flatter PSD, while for higher energies a steeper PSD is found. 
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An important observational test of the model of Comptonization of centrally injected 
soft photons into an inhomogeneous hot medium would be a reliable determination of the 
photon-energy dependence of the PSD power-law slope at high Fourier frequencies. 

In the case of Cyg X-1, Nolan et al. (1981) did not find any significant difference in the 
PSD between the energy channels 12 — 30 keV and 30 — 140 keV, while Brinkmann et al. 
(1974) and Miyamoto & Kitamoto (1989) even found evidence for a flattening of the PSD 
with increasing photon energy, in striking contrast to the prediction of the Comptonization 
model in both geometries investigated here. 

For GX 339-4, no photon-energy dependence of the PSD power-law slope was detected 
( IMotch et al. 1983| , [Maejima et al. 1984|) . 



In the case of an external soft photon source, the early-time slope /i of the single-flare 
light curves also determines the hard time lag (cf. Eqs. [8] and [9]) at high Fourier 
frequencies. The predicted amount of the resulting time lags are comparable for both 
scenarios. The deviation from a power-law behavior of the light curves in the case of an 
central soft photon source, however, leads to a systematic deviation from linearity in the 
dependence of the hard time lag on Fourier period. With the parameters used in Hua et 
al. (1997 b) and in the energy range considered in Miyamoto et al. (1988), this curve is 
well described by a linear dependence, but for higher photon energy channels, a significant 
fiattening of this dependence should be observable if Comptonization of centrally injected 
soft photons is responsible for the rapidly varying hard X-rays. In contrast, such a deviation 
from strict linearity should not be observed in the case of an external soft photon source. 

Very recently, Crary et al. (1997) indeed found a dependence At oc P^'^ for the time 
lag of the 50 — 100 keV photons with respect to the 20 — 50 keV signal from Cyg X-1. 
This seems to be another hint toward the central soft photon injection scenario. Also the 
observed approximate dependence of the Atki oc ln(i?2/-Ei) ([Miyamoto et al. 1988|) is more 
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consistent with this geometry (see Fig. 10 and recall that approximately At oc Act for 
central injection). 

The available data on the phase lags in GX 339-4 and GS 2023+338 (Miyamoto et 
al. 1992) do not show an obvious trend and definitely need to be supplemented by more 
sensitive measurements comparing a larger number of energy channels in order to test the 
Comptonization models investigated in this paper. 

In the discussion of Section 3, we have seen that the break frequencies defining the 
transition between a flat PSD and a steep power-law as well as the transition between the 
linear dependence of the hard time lag on Fourier period to a constant time lag are both 

expected around 

{c/(2nTT R) for central injection 
(12) 
c/ (4 TT R) for external injection 

if Tr >, 1. Since for most Galactic black-hole candidates a Thomson depth of 1 - a few is 
appropriate to model the hard X-ray spectrum via Comptonization of soft photons, Eq. 

(12) yields roughly the same size estimate for both injection geometries, 

{5 • lO^'^ (i^fc-iTr)"^ cm for central injection 
2.5 ■ 10 z/f, _i cm for external injection 

where I'b,-! is the break frequency in units of 0.1 Hz. Since the observed flattening 

frequencies of the PSDs of GBHCs are of the order i/fiat ~ 0.1 Hz and the measurements 
of Miyamoto et al. (1988, 1993) yield a lower limit on the turnover Fourier period of 
the hard time lag in the same frequency range, Eq. (16) reveals a major problem of the 
Comptonization scenario as an explanation for the rapid aperiodic variability of Galactic 
black-hole candidates in general. If indeed the spectral breaks in the PSD and the hard 
time lags are the result of inverse Compton scattering as analyzed in the previous sections, 
then the estimated size of the Comptonizing region is much larger than the expected extent 
of an accretion-disk corona. This problem is obviously even more severe for the case of 
external photons impinging on a hot inner-disk torus. 
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The problem resulting from the break in the power spectrum could be overcome 
assuming a soft radiation flare of much longer intrinsic duration than the light crossing 
time of the inner disk region. In this case, the turnover frequency of the PSD would yield 
an estimate for the typical intrinsic duration of the flare, lyjia^^ ~ 1 — 10 s, which, e. g., is 
of the order of the radial drift time scale in the inner accretion disk. This would generally 
produce an additional break in the PSD at the turnover frequency of the single-shot power 
spectrum. A second break or a smooth steepening of the PSD around 10 Hz has indeed 
been observed in the low states of several GBHCs. 

However, we expect that the assumption of a broad flare as opposed to a 5-shot in 
time does not shift the turnover of the hard-time-lag versus Fourier-period curves if it is 
indeed due to inverse Compton scattering. Therefore, the Fourier-period dependence of 
the hard time lag is still a puzzle in the framework of Comptonization models. As one 
possible solution to this problem we suggest an intrinsic hardening of the soft input photon 
spectrum during the evolution of an extended flare, which could be caused by a hot spot in 
the accretion disk heating up as it spirals inward. This idea will be worked out in detail in 
a forthcoming paper. 



5. Summary and conclusions 

We presented a systematic study on the predicted Fourier power spectra and hard 
time lags resulting from two fundamentally different geometries in Comptonization models 
for the hard X-ray emission of Galactic black-hole candidates: a flaring soft photon source 
located in the center of the Comptonizing region (e. g., a hot corona) and a soft photon 
source located outside a hot inner disk region. The different scenarios yield different 
predictions for the Fourier power spectra and fundamental differences in the photon-energy 
dependence of the Fourier PSD slope and of the slope describing the hard-time-lag versus 
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Fourier-period relation. 

We found that Comptonization of flares of external soft radiation generally leads to 
a weak photon-energy dependence of the power spectrum and the hard time lag, with the 
power spectrum being marginally consistent with the PSD measured for GX 339-1, but 
significantly steeper than for Cyg X-1. Comptonization of centrally injected soft photon 
flares, in turn, leads to a signiflcant photon-energy dependence in the PSD power-law slope 
at high Fourier frequencies and photon-energy dependent deviations from linearity of the 
hard time lags as a function of Fourier period. With a density gradient n{r) oc r^^, the 
simulated PSD for medium-energy X-rays is consistent with the slope observed in Cyg X-1. 

The turnover in the Fourier-frequency dependent hard time lag curve poses a severe 
challenge to Comptonization models for the rapid aperiodic variability in GBHCs in general 
due to the resulting large size estimates. 

Detailed measurements of the energy-dependent properties of the rapid aperiodic 
GBHC X-ray variabihty are strongly encouraged. Using the diagnostic tools developed in 
this paper, they can shed light on the question of the geometry of black-hole accretion disks 
and the location of the Comptonizing region and may serve to rule out a specific geometry 
or even the Comptonization models as a mechanism involved in the rapid aperiodic X-ray 
variabihty of Galactic black-hole candidates in general. 

This work is partially supported by NASA grant NAG 5-4055. 
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Fig. 1. — Eucrgy-dcpcndcnt light curves for Comptonization of an instantaneous soft 
radiation flare at t = in the center of the Comptonizing region with uniform density 
(p = 0) and Thomson depth tt = 1, i? = 3 ■ 10^° cm, kTe = 100 keV 
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Fig. 2. — Energy-dependent light curves for Comptonization of an instantaneous soft 
radiation flare at t — outside a homogeneous Comptonizing region (p — 0) with tt — 3, 
i? = 3 ■ 10^° cm, kTe = 100 keV 
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Fig. 3. — Energy-dependent light curves for Comptonization of an instantaneous soft 
radiation flare at t = outside a homogeneous Comptonizing region with density gradient 
p = 1 with tt = 2.1, = 3 ■ 10^° cm, kT^ = 100 keV 
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Fig. 4. — Energy-dependent power spectrum densities for external soft photon injection, 
p — 0, parameters see Fig. 2 (thick curves: Fourier transform of simulated light curves; thin 
curves: Fourier transform of analytical representation) 
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Fig. 5. — Energy-dependent power spectrum densities for external soft photon injection, 
p — 0, parameters see Fig. 3 (thick curves: Fourier transform of simulated light curves; thin 
curves: Fourier transform of analytical representation) 
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Fig. 6. — Fourier-frequency-dependent hard time lags for internal soft photon injection, 
p — 0; parameters see Fig. 1 
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Fig. 7. — Fourier-frequency-dependent hard time lags for internal soft photon injection, 
p = 1, rr = 0.7, i? = 3 ■ 10^° cm, kT^ = 100 keV 
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Fig. 8. — Fourier-frcqucncy-dependent hard time lags for external soft photon injection, 
p = 0; parameters see Fig. 2. The irregular structures around 0.5 - 1 s are due to numerical 
noise 
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Fig. 9. — Fouricr-frcqucncy-dcpendent hard time lags for external soft photon injection, 
p — 1; parameters see Fig. 3. The irregular structures around 0.5 - 1 s are due to numerical 
noise 
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Fig. 10. — Early-time light curve slopes /i {— a or k, respectively; cf. Sect. 2) as a function 
of the mean photon energy in the different energy bins for the cases p = 0, 1, and 3/2, 
for internal and external soft photon injection. The curves are labeled by the slope density 
gradient slope p and the location of the soft photon source. 
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Table 1. Lightcurve fit parameters for 
Comptonization of a flare of soft radiation in the 
center of the scattering cloud. For definition of 
the parameters see Eq. (1). 





energy channel [keV] 


a 


3 Is] 

y L J 





0.10 - 0.22 


0.65 


0.67 





0.51 - 1.2 


0.75 


0.76 





2.6 - 6.0 


1.55 


0.65 





14 - 31 


2.11 


0.85 





31 - 70 


2.30 


1.05 


1 


0.10 - 0.22 


0.0 


0.77 


1 


0.46 - 1.0 


0.0 


0.71 


1 


2.2-4.6 


0.2 


0.69 


1 


10 - 22 


0.5 


0.63 


1 


46 - 100 


0.6 


0.63 


1.5 


0.10 - 0.22 


- 0.75 


0.90 


1.5 


0.51 - 1.2 


- 0.65 


0.90 


1.5 


2.6 - 6.0 


- 0.80 


0.90 


1.5 


14 - 31 


- 0.80 


0.90 


1.5 


31 - 70 


- 0.80 


0.80 
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Table 2. Light curve fit parameters for Comptonization of a flare 
of soft radiation impinging from outside the scattering cloud. For 
definition of the parameters see Eq. (3). 



p 


energy channel [keV] 


a 


/3[s] 


K 


A/B 





0.10 - 0.22 


1.00 


0.48 


n. det. 


oo 





0.51 - 1.2 


1.00 


0.54 


n. det. 


oo 





2.6 - 6.0 


1.02 


1.55 


n. det. 


oo 





14 - 31 


2.28 


1.30 


1.01 


25.0 





31 - 70 


2.90 


2.00 


1.25 


44.8 


1 


0.10 - 0.22 


8.00 


0.21 


1.00 


9.09 


1 


0.51 - 1.2 


8.00 


0.21 


1.00 


13.5 


1 


2.6 - 6.0 


8.00 


0.25 


1.00 


33.3 


1 


14 - 31 


8.50 


0.29 


1.00 


60.0 


1 


31 - 70 


9.00 


0.32 


2.00 


20.0 


1.5 


0.10-0.22 


4.00 


0.25 


1.00 


4.29 


1.5 


0.51 - 1.2 


4.00 


0.28 


1.00 


11.3 


1.5 


2.6 - 6.0 


8.00 


0.27 


1.02 


13.0 


1.5 


14-31 


8.00 


0.28 


1.00 


14.3 


1.5 


31 - 70 


n. det. 


n. det. 


n. det. 


n. det. 



